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PRINCIPAL  NOTATION 


A Van  Driest  damping  parameter,  see  (2.20b) 

cfg  local  skin  friction  coefficient  in  streamwise  direction 

f transformed  vector  potential  for  ip 

g transformed  vector  potential  for  « 

h^,  h2  metric  coefficients 

h.  net  spacing  in  n-direction 

metric  coefficient,  see  (2.15) 

H boundary- layer  shape  factor,  6*/e^^ 

net  spacing  in  x-di recti on 
geodesic  curvatures,  see  (2.5) 

L mixing  length,  see  (2.20a) 

p static  pressure 

Pg  magnitude  of  velocity  vector  at  the  boundary-layer  edge 

Pj.  Pp  velocity  components  in  boundary-layer  parallel  and  normal, 
respectively,  to  external  streamline 

Rj  Reynolds  number,  u^s^/v 

Rg  Reynolds  number,  Pg®^^/v 

s arc  length  along  coordinate  line 

u,v,w  velocity  components  in  the  x^ , y,  x^  directions 

u^  friction  velocity,  see  (2.20c) 

u^  freestream  velocity 

u^^  reference  velocity 

x^  curvilinear  surface  coordinate 

X2  curvilinear  surface  coordinate  orthogonal  to  the  x-coordinate 

y,  Xj  distance  normal  to  the  body  surface 


x.y.z 


Cartesian  coordinates 


-pu  V , Reynolds  stresses 
-pu'w' 

6 boundary- layer  thickness 

00 

6*  displacement  thickness,  f 0 - 

0 ^ ® 

Gp,  eddy  viscosity 

ejJi  dimensionless  eddy  viscosity, 

n similarity  variable  for  y,  see  (2.22) 


momentum  thickness,  f q./qp(l  — q,./q_)dy 


dynamic  viscosity 
kinematic  viscosity 
density 
shear  stress 

two-component  vector  potential,  see  (2.23a) 


Subscripts 


boundary-layer  edge 
streamwise  direction 


Primes  denote  differentiation  with  respect  to  n 


1.0  INTRODUCTION 


The  work  reported  here  describes  the  progress  made  during  the  past  year 
towards  the  development  of  a general  method  for  computing  three-dimensional 
incompressible  laminar  and  turbulent  boundary  layers  on  ship  hulls.  The 
method  employs  an  implicit  two-point  finite-difference  method  developed  by 
Keller  and  Cebeci^^^  and  an  algebraic  eddy-viscosity  formulation  developed  by 
Cebeci^^^.  During  the'past  year  our  efforts  concentrated  on  the  choice  of  an 
appropriate  coordinate  system;  the  calculation  of  the  geometric  parameters  of 
this  coordinate  system;  the  numerical  solution  of  the  governing  equations  for 
the  orthogonal  curvilinear  system;  and  obtaining  preliminary  results  for 
simple  ship  forms. 

In  Section  2 we  discuss  the  basic  equations  and  the  coordinate  system. 

The  governing  equations  which, in  Section  2.1,  are  written  for  three-dimensional 
incompressible  laminar  and  turbulent  flows  in  an  orthogonal  curvilinear  coord- 
inate system  require  initial  conditions  on  two  intersecting  lines.  For  this 
reason,  in  Section  2.2,  we  discuss  the  formulation  of  the  governing  equations 
to  generate  these  initial  conditions. 

In  principle,  the  calculation  of  boundary-layer  development  on  ship  hulls 
can  be  achieved  in  a number  of  coordinate  systems  applicable  to  arbitrary 
three-dimensional  bodies.  The  presence  of  the  free  surface  suggests  the  use 
of  a streamline  coordinate  system  because  the  flow  boundaries  and  the  coord- 
inate lines  in  that  system  coincide.  In  our  study,  however,  we  abandon  this 

system  in  favor  of  a body-oriented  one  (see  Section  2.2)  suggested  by  Miloh 
(31 

and  Patel'  which  avoids  the  complexity  of  calculating  the  streamline 
coordinates  for  each  Froude  number.  Consequently,  the  coordinate  system  and 
its  geometrical  parameters  such  as  the  metric  coefficients  and  the  geodesic 
curvatures  are  computed  once  and  for  all  for  a specified  ship  hull.  If  the 
ship  hull  is  described  in  analytic  form,  then  these  geometric  parameters  can 
be  computed  exactly.  For  arbitrary  ship  hulls,  they  must  be  computed  numer- 
ically. We  discuss  this  aspect  of  the  problem  in  Appendix  A. 

For  turbulent  flows,  the  governing  equations  require  closure  assumptions 

for  the  Reynolds  stresses.  Here  we  use  the  eddy-viscosity  concept  and  model 

the  Reynolds  shear  stress  terms  by  an  algebraic  eddy-viscosity  formulation 
(21 

developed  by  Cebeci ' . This  phase  of  the  problem  is  discussed  in  Section  2.4. 
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When  physical  coordinates  are  used,  solutions  of  the  governing  boundary- 
layer  equations  are  sensitive  to  the  net  spacings  in  the  streamwise  direction 
(x)  and  the  crosswise  direction  (z),  and  require  a large  nuWber  of  x-  and 
z-stations.  To  remove  this  sensitivity  and  to  generate  the  initial  conditions, 
we  express  the  governing  equations  in  transformed  coordinates.  For  this 
reason,  in  Section  2.5,  we  discuss  a convenient  transformation  and  express  the 
boundary- layer  equations  in  terms  of  transformed  variables. 

In  Section  3 we  describe  the  numerical  method  used  to  solve  the  basic 
equations  for  both  laminar  and  turbulent  flows,  and  in  Section  4 we  present 
preliminary  results  for  a double-elliptic  ship  model.  These  sections  are 
followed  by  a discussion  section  (Section  5)  where  we  outline  the  work  planned 
for  the  next  phase  of  the  studies.  Finally,  in  Appendix  B we  present  a user's 
manual  for  the  computer  program. 


2.0  BASIC  EQUATIONS  AND  THE  COORDINATE  SYSTEM 


2.1  Boundary- Layer  Equations 

The  governinq  boundary-layer  equations  for  three-dimensional  incompressible 
laminar  and  turbulent  flows  in  a curvilinear  orthogonal  coordinate  system  are 
given  by: 


Continuity  Equation 


3x^  ^^2“^  3x2  ^^1’^^  3x3  ° 


(2.1) 


X, -Momentum  Equation 


3x^  hg  3x2  3X3  7 1 h^  3x-|  \p/  3X3  \ 3X3  I 


(2.2) 


X2-Momentum  Equation 

U_3W_+W_3W_+v^+i(u2_.  _1__3_/2.\._3_  / 

h^  3x^  h2  3X2  ^^2  ^ ' ^2  ^^2  ^ ^^3  ) 


Here  h^  and  h2  are  the  metric  coefficients  and  they  are,  in  general, 
functions  of  x^  and  X2,  that  is. 


(2.3) 


h'j  = h^  (x-j  ,X2) 


^2  ~ ^2 ^^1  *^2 ^ 


(2.4) 


The  parameters  K-|  and  K2  are  known  as  the  geodesic  curvatures  of  the  curves 
X2  = constant  and  x^  = constant,  respectively.  They  are  defined  as 


1 ^^1 


1 ^^2 


K = _ ' L 1C.  = - £. 

1 ^1^2  ^^2  h^  h2  3x^ 

At  the  edge  of  the  boundary  layer,  (2.2)  and  (2.3)  reduce  to 


(2.5) 


^ ^ ♦ '=2«e  - Ve  = - R7  5^7  (^) 

u 3W  W 3W„  p „ 

h7  ^ * hj  54  * Ve  - Ve  = - % 0 


(2.6a) 


(2.6b) 


The  boundary  conditions  for  equations  (2.1)  to  (2.3)  are: 


(2.7a) 

(2.7b) 


Xj-O  u,v,w»0 

X3  - 6 u » Ug(x.|,X2)  w » Wg(x^,X2) 

2.2  Initial  Conditions 

The  solution  of  the  system  given  by  (2.1)  to  (2.3)  subject  to  (2.7) 
requires  closure  assumptions  for  the  Reynolds  stresses,  -pu'v'  and  -pv'w' . 
They  also  require  Initial  conditions  on  two  intersecting  planes.  In  general, 
the  construction  of  these  initial  conditions  for  three-dimensional  flows  on 
arbitrary  bodies  such  as  ship  hulls  is  rather  difficult  due  to  a variety  of  bow 
shapes  which  may  be  quite  complicated.  For  this  reason,  it  is  necessary  to  make 
some  assumptions  for  them  in  order  to  start  the  calculations. 

In  our  study  we  choose  the  inviscid  dividing  streamline  on  which  3p/ax2  * 0 
to  be  one  of  the  initial  data  line  (see  Figure  1).  In  the  case  of  a rectilinear 


Figure  1.  A sketch  of  the  initial  data  lines 
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motion  of  a ship,  this  streamline  runs  along  the  plane  of  symmetry.  We  also 
select  this  line  as  the  X2  = 0 line  for  the  coordinate  system.  Because  of  the 
symmetry  conditions  w and  are  zero  on  this  line  causing  (2.3)  to  become 
singular.  However,  differentiation  with  respect  to  X2  yields  a nonsingular 
equation.  After  performing  the  necessary  differentiation  for  the  X2-momentum 
equation  and  taking  advantage  of  appropriate  symmetry  conditions  In  (2.3)  and 
In  the  other  equations,  we  can  write  the  so-called  streamwise  attachment- line 
equations  as: 

Continuity  Equation: 


li,  ("2“)  * syz  * IxjC’iV)  ■ 0 


(2.8) 


-Momentum  Equation: 


u 9u  , au 
h^  ax*|  ^3 


X2 -Momentum  Equation 


u ^»*z  . . .2 


aw. 


K7  5);7  * “z  * ' - '^“'•z  * H 


2 _ ^e  ^'^2^  . . 


aK 


29X2 


U = + w'  - KoU^W_  + 


1 ..2 


ze 


’^‘^e'^ze  ho3Xo  ‘^e 


20^2 


. ax 


3 L 


aw,  


(2.10) 


where  w^  = l/h2  3w/ax2  and  (v'w' )^  = l/h2  3/3x2  (v'w' ).  These  equations  are 
subject  to  the  following  boundary  conditions: 

X2*=0  u»v  = W2  = 0 (2.11a) 

X3  = 6 “ “ “e*  '^z  “ '^ze  (2.11b) 


The  other  Initial  data  line  should  be  selected  near  the  bow  of  the  ship 
along  the  line  perpendicular  to  the  = const,  coordinate.  However,  because 
of  the  variety  of  bow  shapes  possible,  some  of  which  may  be  quite  complicated. 
It  is  necessary  to  make  further  simplifying  approximations  to  the  flow  on  the 
line.  In  our  study  we  assume  the  flow  along  the  Initial  line  to  be  two- 
dimensional  without  pressure  gradient,  and  solve  the  governing  two-dimensional 
equations  for  a flat  plate. 
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2.3  Coordinate  System 

Since,  in  general,  a ship  form  is  a rather  complex  nondevelopable  surface, 
a Cartesian  coordinate  system  is  not  suitable  for  boundary- layer  calculations. 
The  streamline  coordinate  system  is  a very  attractive  one.  The  determination 
of  the  streamlines,  the  orthogonal  lines  and  the  associated  geometrical  param- 
eters require  considerable  work.  In  addition,  they  are  dependent  on  the  Froude 
nuidber  and  the  Reynolds  nunter  if  the  displacement  effect  is  taken  into  account. 
Consequently,  this  coordinate  system,  in  addition  to  being  hard  to  compute, 
becomes  also  uneconomical  to  use  when  the  effects  of  the  Froude  number  and  the 
Reynolds  number  are  to  be  systematically  investigated. 

In  selecting  a coordinate  system  for  the  boundary-layer  calculations,  it  is 
desirable  to  choose  one  such  that  it  is  calculated  only  once.  Although  a gen- 
eral nonorthogonal  coordinate  system  has  been  developed  with  success  by  Cebeci, 
Kaups  and  Ramsey the  numerical  work  for  evaluating  the  geometrical  param- 
eters is  very  involved.  For  this  reason,  in  the  present  study,  we  have  decided 
to  abandon  this  coordinate  system  in  favor  of  the  one  suggested  by  Miloh  and 
Patel This  coordinate  system  depends  only  on  the  geometry  of  the  body  and  is 
calculated  once  and  for  all.  Therefore,  the  effects  of  the  draft,  the  trim,  the 
Froude  nurrtber,  and  the  Reynolds  nuiflber  on  ship  boundary-layer  development  can 
be  systematically  studied.  In  this  coordinate  system,  one  of  the  coordinates 
is  taken  as  the  lines  of  * constant.  The  other  coordinate  X2(x,z) 
constant,  which  is  orthogonal  to  x^  = constant  lines  everywhere  on  the  ship 
hull,  can  be  shown  to  be  the  solution  of  the  differential  equation 


dz 

ax 


^x^z 

1 +7? 


(2.12) 


Here  y = f(x,z)  defines  the  ship  hull,  and  (x,y,z)  denote  the  Cartesian 
coordinates.  The  geometrical  parameters  are  calculated  from  the  following 
formulas 


0 ^ 

^2  3X2/3Z 


(2.13a) 


(2.13b) 
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''i  =- 


(2.13c) 


%^zz 

(1  + (1  + f^  + f^)^ 


(2.13d) 


It  may  be  noted  that  K-j  and  K2  can  also  be  obtained  from  ( 2.5)-  This 
provides  a check  on  the  expressions  given  above. 

In  the  boundary-layer  equations,  h2  always  occurs  In  the  combination 
l/h2  3/3x2.  Since  we  have  chosen  x and  X2  as  Independent  variables,  we 
have 


where 


1_  _3_  . ]_  3Z 3_  , 1 _3_ 

^2  ®^^2  ^2  ^^^2  ^^z  ^ 


4 3Z 


(2.14) 


(2.15) 


From  (2.13b)  it  follows  that 


h^  = (1  + f^)’^^ 


(2.16) 


We  see  that  once  the  curvatures  have  been  found.  It  Is  no  longer  necessary  to 
work  with  (x,  X2)  coordinates;  instead  we  can  work  with  (x,z)  coordinates 
layers  on  ship  hulls  and  solve  the  governing  equations  In  that  system.  The 
use  of  this  coordinate  system  offers  the  major  computational  advantage  that 
the  pressure  distribution  can  be  prescribed  much  more  readily  In  terms  of 
X and  z than  In  terms  of  x and  X2. 

In  the  boundary- layer  calculations,  we  need  the  velocity  components  along 

the  surface  coordinates.  For  the  coordinate  system  adopted.  It  can  be  shown 

(5) 

that  the  velocity  components  In  the  two  systems  are  related  by'  ' 


(1  + fpu  ^ y - f/^W 

[(1  + f2)(l 


(2.17a) 
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Here  (U,  V.  VI)  are  the  velocity  components  In  the  Cartesian  coordinates, 
and  (Ug,  Mg)  In  the  adopted  coordinates. 

2.4  Turbulence  Hodel 

For  turbulent  flows.  It  Is  necessary  to  make  closure  assumptions  for  the 
Reynolds  stresses,  -pu'v'  and  -pv‘w‘ . In  our  study,  we  satisfy  this 
requirement  by  using  the  eddy-viscosity  concept  and  relate  the  Reynolds  stresses 
to  the  mean  velocity  profiles  by 


-u'V 


^m  3y 


^m  3y 


(2.18) 


Here  we  have  let  Xg  = y.  Vie  use  the  eddy-viscosity  formulation  of  ref.  2 and 
define  by  two  separate  formulas.  In  the  Inner  region,  Is  defined  as 


where 

L = 0.4y  [1  - exp(-y/A)] 


A = 26vu'^ 

T 


In  the  outer  region,  Ej^  Is  defined  by  the  following  formula: 


(2.19) 


(2.20a) 

(2.20b) 

(2.20c) 


(2.21) 


The  Inner  and  outer  regions  are  established  by  the  continuity  of  the  eddy 
viscosity  formula. 


2.5  Transformation  of  the  Basic  Equations 


The  boundary- layer  equations  can  be  solved  either  in  physical  coordinates 
or  in  transformed  coordinates.  Each  coordinate  system  has  its  own  advantage. 
In  three-dimensional  flows,  the  computer  time  and  storage  required  Is  an 
important  factor.  The  transformed  coordinates  are  then  favored  because  the 
coordinates  allow  larger  steps  to  be  taken  in  the  streamwise  and  spanwise 
directions. 


We  define  the  transformed  coordinates  by  (x-  = y) 

1/2  X, 


xi  = x^. 


X2  = 


dn 


■ft)" 


h^dx^ 


(2.22) 


and  introduce  a two-component  vector  potential  such  that 


wh  = ^ 
'^"1  ay 


wk  h = — / 

''"l"2  ^ax^  axg/ 


(2.23a) 

(2.23b) 

(2.24a) 

(2.24b) 


where  (|>  and  41  are  defined  as 

4)  - (vUgS^j ) ^ h2^(Xi  ,X2 ,n) 

<l>  " (vUgSi ) ^ *^ref^*^e 

and  u^g|r  is  some  reference  velocity. 

Using  these  transformations  and  the  relations  given  by  (2.6)  and  (2.12), 
after  some  rearranging,  we  can  write  the  x -momentum  and  z-momentum  equations 
for  the  general  case  as: 

x-Momentum 


(bf)'  + m^ff  -m2(f)^  - V‘9'  + ->"8^9')^  + 


m 


11 


= m 
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z-Momentum 


(bq")'  + m^fg"  -tn3(g')^  + nigqg''  -mg(f  + m^2 

= ■ ""  fj^)  * "7  (s'  “ 0"  1^) 


(2-26) 


and  their  boundary  conditions  as 


n = 0,  f = f = q = g'  = 0,  n = n„.  f = 1,  g'  = 

(2.27) 


Here  primes  denote  differentiation  with  respect  to  n,  and 
f = u/Ug,  g'  = w/u^^,  b = 1 + 4’  4 = 

The  coefficients  m^  to  m^2  given  by 


2 ■ Ug  3S^  ’ 


m3  = 0 


"’4  “ "^1*^2’ 


, _ Vef  h !!!e  , , 

^ - Ug  Ugh2  3X2  ^ ^ 


-Vf  3 J ^ .!Vef!i 

- Ug  [2^3x2  Ug  h23x2/  ”'7  Ug  h2 


■ ( If)  • “a  = 'I'^l  (u^J  • 


"’ll  " ^1 


L_  !!!e  + !e  /!_  _?!!e_  ^ ^ ^-K^l 
“e  "'e  W 2 Ug  N/J 


m - 

■"lo  • R7 


_!!§_  f!i  !!e  WM  !i 
“ref  K \“e  / “e 


3W  W 

fusxT  ^ uT 

c c e 


(2.28) 


To  transform  the  streamwise  attachment-line  flow  equations  and  their 
boundary  conditions,  (2.8)  to  (2.11),  we  define: 


""2 ' 1]^  • 


1 w - ^ 
’l“z  3y 
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V 


(2.30b) 


' - * '‘2*) 

with  t|)  and  ({>  still  given  by  (2.24).  With  these  variables,  the  streaitwise 
attachment-line  equations  in  the  transformed  coordinates  can  be  written  as 

(bf-)'  = ”10  (^' 

(bg")  + m^fg"  -m^f'g'  -m3(q')^  + mggg"  -mg(f  )^  + m^2 

The  boundary  conditions  and  the  coefficients  m-j  to  m^2  same  as 

those  in  (2.27)  and  in  (2.29)  except  now 


n -!iVf 

"3-h2  u^ 


m 


n = n. 


aig  — nij , mQ  — 


g'  = w /u 


ze'  “ref 

3K, 


(2.33a) 


_ . ^1  '“e 

'9  ■ u^g^  "1  h23X2  ’ "’ll  ■ Ug  3s^ 


(2.33b) 


12  " u 


ref 


!i!r!se,,  (V)  _SK 


The  governing  equations  for  the  initial  conditions  normal  to  the  stream- 
wise  attachment-line  flow  are  given  by  the  two-dimensional  flat-plate  equations, 
which  in  terms  of  transformed  variables  are 

(bg")'*^fg"  = m,(g'|^-g"|a-j 

For  these  equations,  the  boundary  conditions  are  given  by  (2.27). 

In  terms  of  transformed  variables,  the  algebraic  eddy-viscosity  formulas 
as  given  by  (2.13)  to  (2.15)  become 


(2.34) 

(2.35) 
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3.0  NUMERICAL  METHOD 
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We  use  the  Box  method  to  solve  the  governing  equations.  This  is  a 
two-point  finite-difference  method  developed  by  Keller  and  Cebeci.  So  far 
this  method  has  been  applied  to  two-dimensional  flows  as  well  as  three- 
dimensional  flows  and  has  been  found  to  be  efficient  and  accurate.  Descrip- 
tions of  this  method  have  been  presented  in  a series  of  papers  and  reports, 
and  a detailed  presentation  is  contained  in  a forthcoming  book  by  Cebeci  and 
Bradshaw^®\  For  completeness  we  shall  present  a brief  description  of  this 
method  in  application  to  our  equations.  At  first  we  shall  discuss  the  solu- 
tion of  the  streamwise  attachment-line  equations,  (2.31)  and  (2.32),  and  then 
the  solution  of  the  full  three-dimensional  flow  equations  as  given  by  (2.25) 
and  (2.26). 


3.1  Difference  Equations  for  the  Streamwise  Attachment-Line  Equations 

According  to  the  Box  method,  we  first  reduce  the  equations  (2.31),  (2.32), 
(2.27)  and  (2.33a)  into  a system  of  six  first-order  equations  by  introducing 
new  dependent  variables  u(x^,  n).  v(x^,  x^,  n).  w(x^  X2»  n)  and 

t(x^ , X2,  n).  Equations  (2.31)  and  (2.32)  then  can  be  written  as 

f = u 


u'  = V 


g'  = w 


w'  = t 


(bv)'  * m,f.  - n,/  * mgg»  . m„  ■ n>,(,  (u 
(bt)'  ♦ - m,uw  - mjW^  t m^gt  - ♦ ",2  = ”,0  ("  ‘ 


The  boundary  conditions  (2.27)  and  (2.33a)  become 
0*0  f = q = u = w = 0 

n = n u = 1 


= «ze/“ref 


(3.1a) 
(3.1b) 
(3.1c) 
(3. Id) 

(3.1e) 


(3.1f) 


(3.2a) 

(3.2b) 


L 


We  next  consider  the  net  rectangle  shown  in  Figure  2 and  denote  the  net 
points  by 


(x,)o  = 0 

(*1  >n  ' ("l )„-l  * "n 

n » 1 

o 

II 

o 

"J  ’ "M  * "j 

J « 1 

n 


We  approximate  the  quantities  (f,  u,  v,  g,  w,  t)  at  points  nj] 

of  the  net  by  functions  denoted  by  (fj.  Uj,  Vj,  g”,  Wj,  tj).  We  also 
employ  the  notation  for  points  and  quantities  midway  between  net  points  and 
for  any^  net  function  Sj : 

^^l^n-1/2  “ l^^^l^n  * ^^l^n-1^  '’j-1/2  “ T ^j-1^ 

(3.3) 

n-1/2  ^1  / n 1 /c'’  a \ 

r ' *j-l/2  7 ^j-1^ 

The  difference  equations  which  are  to  approximate  (3.1)  are  formulated  by 
considering  one  mesh  rectangle  as  in  Figure  2.  We  approximate  (3.1a,  b,  c,  d) 
using  centered  difference  quotients  and  average  them,  about  the  midpoint 
[(xi)nt  of  the  segment 

~ ^j-1^  “ “j-1/2 
“ “j-1^  " ''j-1/2 
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^ - '^j-1/2 


(3.4c) 


h'^  (w"  - w"  ) = t” 
'^j-V  ^j-1 


(3.4d) 


Similarly,  (3.1e,f)  are  approximated  by  centerino  them  about  the  midpoint 
(x^)n_i/2.  rectangle  P-|P2^3P4-  This  gives: 


Kj'[(b.)5  - (bv>;.,]  * (mj  . - (mj  * .„)(u^)3_,/2 


+ ">6(gv)j.y2  “n'^''j-1/2^j-l/2  ^j-1/2''j-l/2J 


•^j-1/2  "'ll 


(3.4e) 


hj^[(bt)J  (bt)J_i]  + (tn!J  + «n)(Tt)J_i/2  (mj  + “n^^"'^^j-l/2 


- m3(w  )j_^^2  ■'■  "'6^9t)j.i/2  "’9^“  ^j-1/2 


+ a Fw'^'^  u"  - u""^  w"  + t""^ 

^ “n'-'^j-l/2''j-l/2  "j-l/2'^j-l/2  ^j-l/2^j-l/2 


^n-l  -n  T _ T-n-1  „n 
j-l/Z^j-l/ZJ  ■ 'j-1/2  '"12 


(3.4f) 


j-l'/2  “ ^ j-1/2  ^“^^j-1/2^  '"l  j-1/2 


“1  wu..\n-l 


^n-l  I , _n-l  vn-1 


n-l/..2,n-l  X X mO-li 

2 ' ^j-1/2  "^6  ^^''^j-1/2  11  J 


^j-i/2  " ^ j-1/2  ^“‘^^j-1/2^  ^ ""l  '^^^^j-i/2 


(3.5a) 


^n-1 I , _n-l/^. \n-l 


">4  (“w)j.y2  "’3  ^j-1/2  "’e  ^^^^j-1/2 

n-l,..2vn-l  X 
9 ^j-1/2  "'12  ^ 


(3.5b) 
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m. 


° - (*1  'a-l 


n-1/2 
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3.2  Difference  Equations  for  the  Full  Three-Dimensional  Equations 


(3.5c) 


The  difference  equations  for  the  full  three-dimensional  equations,  as 
given  by  (2.25)  and  (2.26), are  again  expressed  in  terms  of  a first-order 
system.  With  the  definitions  given  by  (3.1a)  to  (3. Id),  they  are  written  as 


(bv)'  + m^fv  - m2U  - m^uw  + mggv  - mgW  + m^^ 


(3.6a) 


(bt)'  + m^ft  — m^uw  — m^w  + mggt  — m^u  + 


m 


= m. 


"10  ^3X^)‘^"'7(” 

Their  boundary  conditions,  (2.27),  become: 
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3JV__  . 


3X, 


3X, 


) 


(3.6b) 


n = 0 


n = n. 


f=g=u=w=0 


u-l. 


(3.7a) 

(3.7b) 


The  difference  equations  for  (3.1a)  to  (3. Id)  are  the  same  as  those  given 
by  (3.4a)  to  (3.4d);  they  are  written  for  the  midpoint  [(x^)^^,  '^j-1/2^ 

of  the  net  cube  shown  in  Figure  3;  that  is, 

‘'j-1/2  ’ ''j-1/2’ 

The  difference  equations  which  are  to  approximate  (3.6a ,b)  are  rather 
lengthy.  To  illustrate  the  difference  equations  for  these  two  equations,  we 
consider  the  following  model  equation 


(bv)' 


(3.8) 
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n(j)  (j.n,i-l) 


(j.n.i) 


”T  (j-T,n, 

1 ,4^ 


(j-r.n-M) 


Figure  3.  Net  cube  for  the  difference  equations  for  three-dimensional  flows. 
The  difference  equations  for  (3.8)  are: 

~ ^^j-1^  * ^'"l^j-l/2^^^j-l/2  ^ ^"’ll^j-1/2 


/_  \n-l/2  — n II- 1 , /_  \ii-i/£.  — I 

" '"’lOM-1/2  *^^-1/2  ^"’7M-l/2  '^j-1/2  r. 


n n- 1 


n-1/2  - 


“i  -^-1 


Here,  for  example. 


%=i(r 

^ ‘'j  ^ “j-1  ^ “j-1  ^ 

(m,  Mm,)?-’  Mm, )J;|] 


(3.9) 


(3.10) 


^0  = 0 


^i  = "i-1  " "i 


i = 1,  2,  ...,  I (3.11) 
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3.3  Solution  of  the  Difference  Equations 

The  difference  equations  (3.4)  for  the  streatnwise  attachment-line  flow 
and  the  difference  equations  for  (2.34)  and  (2.35)  toqether  with  the  differ- 
ence equations  for  (3.1)  are  nonlinear  algebraic  equations.  We  use  Newton's  methc 
to  linearize  them  and  then  solve  the  resulting  linear  system  by  the  block- 
elimination  method  discussed  by  Keller  (1974).  A brief  description  of  it  will 
be  given  for  the  streamwise  attachment-line  equations. 

Using  Newton's  method,  the  linearized  difference  equations  for  the  system 
given  by  (3.4)  and  (3.2)  are: 


h. 

6U.  -6Uj-i  -^(6V.  . 6V..^)  = (r^). 

h. 

6g^-6<Jj.l  -/(6W.  + 6W..^)  = (r3). 

h. 


(3.12a) 

(3.12b) 

(3.12c) 

(3.12d) 


,w.  -6Wj_i  + 6t..,)  - (r^). 

(Ci)j4v.  + (^2)j«Vj.i  + (C3)j.6f.  + (^4)j6f._^  + (5s)j6U.  + (C6)j.6U..^ 

+ Uy)^6g.  + (CsV^j.!  = (^)j  (3.12e) 

(6l).6t.  ^ (e2)j«t._i  + (33)j6f.  + (B4).6f..i  + (B5).6w.  ^ (B6)j6w._^ 

+ (B7)j6u.  + (B8)j.4Uj_^  + (3g).6g.  + <3.12f) 

Here  we  have  dropped  the  superscripts  n,  i and  have  defined  (r. ).,  (?. ). 

K J K J 

and  (6|^)j  by 


= 'j-l  • 


- "j-1  -“J  * "j‘,1-1/2 


(3.13) 
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'j-1/2  "’ll 


[(bv)j_i/2  + (m^  + On)(^v)j  1^2  (">2  ■*■  “n^ ^”^^--l/2 

+ "’6^‘’''Vv2  ^ “n  ^''j-l/2^j-l/2“  ^j-l/2''j-l/2^^ 

i-1/2  “"’12  “ ^^'^^^j-1/2  ^"’1  ■’■  “n^^^^^j-1/2  “ ^'"4  ■’■  “n^^“‘^Vl/2 

”"’3(w^)j-l/2  ■’■  "'6^^^Vi/2  ~"’9^“^Vi/2 


* “n<"3-l/2“j 

-1/2 

: “j-l/2'^j-l/2  ^ 

s": 

1 f 

l/2^j- 

-1/2 

“ ^j-l/2^j-l/2^^ 

(3.13) 

(i;,)j  = 

b.  , 

^ 2 ^""l  ^ “n 
J 

“n  ,n-l  1 

2 ^j-1/2  2 *"( 

i’j 

(3.14a) 

<^2>i  = 

b.  , , 

- ^ 2 ^"’1 

<3 

■*■  a 

n 

)f  .^  !a  f^-l 

^ j-1  2 ^j-1/2 

Vj- 

(3.14b) 

f <"’1  * “n’^j  * 

“n 

2 

J-1/2 

(3.14c) 

<«4>J  ' 

7 ("1  * '■„)’'j-1 

a 

+ - 

2 ''j-1/2 

(3.14d) 

-(m2  k ,„)Uj 

(3.14e) 

-(m2  t a„)Uj_, 

(3.14f) 

*‘7^k  ° 

? "'e'-j  > 

"’fi 

(^8>j  = r »j-i 

(3.14q) 

(C,)j. 

(62>j  = (cjjj. 

(B3)j.  = 2 S “n^^j 

t"-'' 

^j-1/2 

(«4>J 

' i *"l  * 

■1  ^ 

“n  .n-1 
2 ^j-1/2’ 

<«5>j 

* " 7 <"’4  * “n>“j  ~ 

°*n  n-1 

"’3’^j  “ “2  ’’j-1/2 

<“6)j 

'-7<"4*“n>“j-l 

“"'3’^j-l  “T”j- 

1 

1/2 

(3.15) 

(Py)j  = 2 K ^ Vj  * 2 ’^j-1/2 

= ~ F K “n^'*j-l  ~ Vj-1  T”j-i/2 

1 nif 

(69)j  = 2Vj*  (B-lO^j  * 2"  ^j-1 
The  boundary  conditions  becone 

5fo  = 69q  = = 0 ♦ BUj  = 6Wj  ' 0 (3.16) 

The  solution  of  the  linear  system  given  by  (3.12)  and  (3.16)  is  obtained 
by  using  the  block  elimination  method.  According  to  this  method,  the  system 
is  written  as 


' S 

Here 


(3.17) 


Cl 

(r,)j 

ro 

6Uj 

Cj 

ill 

s-<? 

(r4>j 

ill 

69j 

6W. 

1 

1 

The  Aj,  Bj,  Cj  in  A denote  6x6  matrices.  The  solution  of  (3.17)  is 
obtained  by  the  procedure  described  in  reference  6. 
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4.0  PRELIMINARY  RESULTS 


In  order  to  check  the  method,  we  have  considered  the  calculation  of 
boundary- layer  development  on  a body  represented  analytically  by 


y = f(x,z)  = B 

’-(i 

1/2 

’-(i 

1/2 


(4.1) 


This  is  a double-elliptic  ship  having  round  edges  except  for  the  shape  corners 
at  X » +L  and  z * +H  as  shown  In  figure  4.  The  potential -flow  solutions 
for  the  body  of  L:H:B  = 1:0.1:0.125  were  obtained  from  the  Douglas-Neumann 


Figure  4.  The  double-elliptic  ship  given  by  (4.1). 


computer  program  for  three-dimensional  flows  (see  figure  5).  To  get  the  solu- 
tions, 120  control  elements  on  the  surface  were  used,  12  along  the  x-directlon 
and  10  along  the  z-dl recti on. 

Before  we  describe  our  boundary- layer  calculations  for  the  ship  given  by 
(4.1),  It  Is  useful  to  discuss  the  pressure  distribution  for  this  body.  As 
can  be  seen  from  figure  5,  the  streamwise  pressure  gradient  Is  Initially  favor- 
able In  the  bow  region  and  then  adverse  up  to  the  midpoint  of  the  body.  This 
is  followed  by  a region  of  favorable  pressure  and  then  by  a sharp  adverse 
pressure  gradient  very  close  to  the  stern.  The  crosswise  pressure  gradient 
varies  In  a more  complex  manner.  Near  the  bow,  the  pressure  decreases  as  z 
Increases.  As  the  flow  moves  downstream,  the  location  of  the  minimum  pressure 
moves  up  and  finally  reaches  the  water  surface.  This  Implies  that  the  cross- 
flow  will  reverse  Its  sign  a number  of  times.  As  a result.  In  some  regions  we 
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may  expect  the  cross-over  type  of  the  crossflow  profiles  (I.e.,  where  the 
direction  of  the  crossflow  changes  across  the  boundary  layer).  This  conclusion 
is  drawn  from  considering  the  pressure  gradients  only.  The  real  situation  is 
more  complicated  because,  in  addition,  there  are  the  upstream  effects  and  the 
curvature  effects  on  the  flow  characteristics. 

The  boundary-layer  computation  starts  from  X'j/L  ■ -0.97.  The  region 
between  this  point  and  the  bow  is  not  considered  for  two  reasons:  first,  the 
pressure  distribution  in  this  region  is  not  available  from  the  potential -flow 
computer  program  and  second,  the  sharp  corner  at  x^/L  = -1.0  and  z/H  ■ 1.0 
implies  an  infinite  curvature  and,  hence,  the  boundary-layer  assumptions  are 
no  longer  valid  in  that  region.  In  our  boundary- layer  calculations,  we  have 
used  32  points  along  the  x^-direction  and  12  points  along  the  Xg-direction. 

In  the  normal  direction,  we  have  taken  approximately  40  points.  The  flow  was 
assumed  to  be  laminar  at  first  and  then  was  specified  to  be  turbulent  at  the 
third  (x^/L)-station,  namely  at  x^/L  » -0.925. 

Some  of  the  computed  results  for  R|^  = 10^  are  shown  in  figures  6 to  10. 
Figures  6,  7,  and  8 show  the  streamwise  variation  of  the  local  skin-friction 
coefficient,  the  shape  factor,  and  the  Reynolds  number  based  on  the  momentum 
thickness  for  z/H  = 0,  0.1,  0.2,  and  0.3,  respectively.  Figures  9 and  10 
show  some  of  the  typical  streamwise  and  crosswise  velocity  profiles  along 
z/H  = 0.3.  In  these  two  figures,  the  profiles  for  z/H  = 0.0  correspond  to 
the  streamwise  attachment  solutions  (z/H  = 0).  Vfe  see  from  figure  3 that  the 
flow  along  the  attachment  line  separates  at  approximately  x/L  - 0.95,  or 
near  the  stern,  where  the  pressure  gradient  is  sharply  adverse.  As  z/H 
increases,  the  point  of  separation  moves  upstream.  Along  the  line,  z/H  - 0.3, 
the  flow  separates  at  about  x/L  = 0.80,  where  the  pressure  is  almost  constant. 
This  unusual  behavior  results  from  the  crossflow  effect.  A better  marching 
procedure  is  required  to  remedy  the  situation  as  is  discussed  in  the  next 
section. 


Figure  6.  Variation  of  skin-friction  coefficient  alonq  the  streamwise  direction  for  the  double-ell iotic  ship. 


Figure  8.  Variation  of  along  the  streamwise  direction  for  the  double-elli 
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5.0  DISCUSSION 


The  present  effort  is  one  step  towards  developing  a general  method  for 
computing  three-dimensional  boundary  layers  on  ship  hulls.  Obviously  further 
studies  are  needed  in  the  following  areas: 

1.  A better  numerical  procedure  for  handling  the  computation  of  flows 
in  which  there  are  regions  of  negative  crossflow.  Our  present  cal- 
culations here  and  in  the  related  problems^ ^ ^ demonstrate  that  in 
regions  of  the  flow  field  where  there  is  a significant  layer  of 
crossflow  opposing  the  marching  direction,  oscillations  can  occur 

in  the  generated  solution.  These  are  caused  by  the  stability  condi- 
tion on  the  difference  equations  which  does  not  permit  the  correct 
domain  of  dependence  of  the  differential  equation  to  be  utilized  at 
the  particular  point  being  computed. 

2.  A better  coordinate  system  than  the  one  used  in  this  study  is  required. 
As  discussed  In  Section  2.3,  a nonorthogonal  system  Is  appropriate  for 
ship  hulls.  Although  the  governinq  equations  are  slightly  more  com- 
plicated than  those  for  an  orthogonal  system,  that  complication  does 
not  cause  any  difficulties  in  their  numerical  solution.  The  big 
advantage  of  the  nonorthogonal  system  is  the  control  it  allows  over 
the  coordinate  net  spacing.  In  orthogonal  systems,  a family  of 
surfaces  which  intersects  the  hull  surface  has  been  selected,  the 
orthogonal s are  consequently  determined  and,  for  some  ship  hulls, 

can  lead  to  coordinate  nets  which  are  unacceptably  sparce  at  some 
areas. 

These  necessary  studies  indicated  above  are  presently  in  progress.  A 
recently  developed  procedure,  which  looks  particularly  promising,  uses  the 
notation  of  domain  of  dependence  and  follows  the  characteristics  of  the 
locally  plane  flow.  In  this  procedure,  the  direction  of  w at  each  nj-grid 
point  is  examined  and  the  difference  equations  formulated  accordingly.  In 
addition,  the  new  procedure  utilizes  a modification  of  the  Box  scheme  similar 
to  the  zig-zag  differencing  scheme  proposed  by  Krause  et  al.^®^  and  success- 

/ Q \ 

fully  used  by  Cebeci'  ’ to  compute  the  unsteady  flow  generated  by  an  impul- 
sively started  circular  cylinder.  The  finite-difference  equations,  with  the 


! 
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zig-zag  differencing  used  to  advance  in  time,  have  been  solved  from  time  t = 0 
until  the  time  the  separation  point  has  moved  from  the  rear  of  the  cylinder  to 
its  final  steady-state  value  of  e = 106“.  The  flow  reversal  has  been  accounted 
for  satisfactorily  and  the  calculations  extended  further  in  time  than  many 
previous  studies. 

An  example  of  a nonorthogonal  coordinate  system  is  given  by  the  intersec- 
tion of  meridian  planes  and  parallel  cuts.  Consider  the  hull  as  given  in  the 
usual  Cartesian  coordinate  system,  that  is,  x along  the  ship  axis,  y and 
z in  the  crossplane  (see  Figure  11).  We  select  x as  one  of  the  coordinates 
and  the  polar  angle  ^ in  the  yz-plane  as  the  other  coordinate.  It  is  only 
necessary  to  define  the  hull  as  a family  of  points  in  the  planes  x = constant. 
Then  the  data  is  spline-fitted  in  each  crossplane  as  y vs  (|i  and  z vs  (t*, 
and  interpolated  for  constant  values  of  <|i.  Then  another  set  of  spline  fits 
in  planes  (j>  = constant  for  y vs  x and  z vs  x completes  the  definition  of 
the  coordinate  system.  The  lines  formed  by  the  intersection  of  the  planes 
X = constant  and  » constant  with  the  hall  surface  constitute  the  nonortho- 
gonal coordinate  net  on  the  surface,  and  the  third  boundary- layer  coordinate 
is  taken  as  the  distance  normal  to  the  surface  in  accordance  with  first-order 
boundary- layer  approximations. 

Since  the  spline-fitting  also  yields  derivatives,  the  metric  coefficients 
and  the  geodesic  curvatures  of  the  coordinate  lines  can  be  calculated  by  using 
the  formulas  given  below. 

The  metric  coefficients: 


The  ;inn1p  hptwppn  thp  rnnv>HiriAtp  11np< 


The  geodesic  curvature  of  the  = const,  line 

[(w)x  (^l  ■ (^)x]  • ^ » f?) 

' <p 

" [(^)x 

The  geodesic  curvature  of  the  x = const,  line 

The  other  curvature  terms  are  calculated  from  the  following  expressions  given 
in  reference 

"12  ' iW  [-  ("1  ^ B7fr)*  (l<2  * 
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APPENDIX  A 


CALCULATION  OF  GEOMETRIC  PARAMETERS 
OF  THE  COORDINATE  SYSTEM 

As  discussed  in  Section  2.3,  the  geometrical  parameters  and  external 
velocity  components  along  the  adopted  coordinates  require  the  quantities 
fx»  f^.  etc.,  on  a ship  hull.  A ship  hull  is  usually  described  in  the  form 
of  a family  of  curves.  For  example,  it  may  be  given  as  the  curves  generated 
by  the  intersection  of  the  hull  with  planes  (x  = constant)  perpendicular  to 
the  free  surface  at  various  stations  along  the  ship  hull.  The  problem  is 
then  to  find,  with  sufficient  accuracy,  the  derivatives  at  an  arbitrary  point 
on  the  hull.  When  the  ship  hull  is  described  analytically,  then  these 
quantities  and  consequently  the  geometric  parameters  of  the  coordinate  system 
can  be  computed  exactly.  For  arbitrary  ship  hulls,  they  must  be  computed 
numerically. 

The  exact  solutions  for  the  three-dimensional  velocity  potential  is 
not  available  even  on  a double  ship  model,  where  the  free  surface  is  treated 
as  a rigid  plane.  It  is  usually  given  experimentally  or  calculated  numerically 
by  the  method  of  Hess  and  Smith.  In  this  method  the  singularities  (e.g.,  source, 
dipole,  vortex,  etc.)  are  distributed  on  the  body  surface,  and  the  governing 
equation  is  transformed  from  Laplace's  equation  in  the  field  into  a Fredholm 
integral  equation  over  the  body  surface.  This  integral  equation  expresses 
the  condition  of  zero  normal  velocity  on  the  surface.  To  solve  the  integral 
equation  numerically  the  body  surface  is  divided  into  a large  number  of  small 
surface  elements,  which  may  be  plane  or  curved.  On  each  element  the  strength 
of  a singularity  is  assumed  to  be  either  constant  or  variable  in  such  a way 
that  the  overall  distribution  may  be  defined  by  a single  unknown  parameter  for 
each  element.  The  normal -velocity  boundary  condition  is  required  to  hold  at 
one  so-called  control  point  on  each  element.  This  yields  a number  of  linear 
equations  equal  to  the  number  of  unknown  parameters  defining  the  singularity 
distribution.  This  set  of  equations  is  solved  by  either  direct  elimination  or 
Gauss-Seidel  iteration. 

The  direct  elimination  is  accomplished  by  the  Douglas-developed  SOLVIT 
routine,  which  is  applicable  to  a large  full  matrix  that  will  not  fit  into 
high-speed  core  storage.  Recently  this  method  has  been  shown  to  be  much 
faster  than  standard  IBM  subroutines.  Once  the  equations  have  been  solved,  the 
resulting  singularity  distribution  is  then  used  to  calculate  values  of  fluid 
velocity  and  pressure  at  the  control  points.  As  can  be  easily  understood  the 
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control  points  of  the  potential -flow  computations  will  not  coincide  with  the 
points  whei^e  the  boundary  layer  Is  to  be  computed.  There  Is,  therefore,  a 
need  to  develop  a numerical  method  to  Interpolate  a function  of  two  Independent 
variables. 

There  are  a number  of  well-known  numerical  techniques  that  can  be  used  for 
Interpolation  and  differentiation  of  a function  of  one,  two,  or  more  variables. 
Here  we  shall  make  use  of  the  surface  splines  for  Interpolating  and  obtaining 
the  first  derivatives  of  a function  of  two  variables  and  the  Fourier  series 
expansion  for  finding  the  second  derivatives. 

The  surface  spline  Is  based  on  the  small  deflection  equation  on  an  Infinite 
elastic  plate  that  deforms  In  bending  only.  The  procedure  Is  to  represent  the 
deflection  due  to  a point  load.  The  deflection  of  the  entire  spline  Is  then 
taken  as  the  sum  of  all  the  point  load  distributions  subject  to  the  boundary 
condition  that  the  surface  should  become  flat  a long  distance  from  the  applied 
loads.  This  results  In  a system  of  linear  equations  which  Is  solved  for  the 
spline  coefficients.  The  final  system  of  equations  Is  given  below.  (For 
details  of  derivation,  see  ref.  10.) 


II 

En 


(A.l) 


(A.  2) 


IV 

E = 


(A.  3) 


IV 

•o  * * vj  *E  ■'ij  ■ "j-  ^ 


(A.  4) 


where 


"ij  = 
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If  the  function  being  Interpolates  has  one  or  two  planes  of  symnetry  or 
antisymmetry,  then  use  can  be  made  of  Images  to  Improve  the  accuracy  of  the 
fit  or  to  reduce  the  number  of  simultaneous  equations.  If  w Is  symmetrical 
about  x^,  the  system  of  equations  becomes 

N 

1-1 
N 

23  "1^1  * “ 

1-1 

N 

*0  * Vj  1"  = Wj  j = 1 ,N  <*•«> 

1=1 

where 

» (X,  *Xj-2Xo)t  (Zj-Z,)2  „ 

Similar  results  can  be  obtained  for  z symmetry  and  for  both  x and  z 
symmetry.  The  system  of  equations  (A.l)  to  (A. 5)  Is  solved  using  the 
above-described  SOLVIT  routine. 

Equation  (A. 4)  can  be  differentiated  to  give  the  first  derivatives  of  w, 
or  the  slopes  of  the  spline 

, N 

9W.  p 

^ = ai  + 2 ^ F^(l  + In  r^j)(xj  - x^) 

. N 

3^-  ag  +2  2^  Fid  + In  r^j)(Zj  -z^) 

1=1 

Though  (A. 4)  can  be  successively  differentiated  to  give  higher-order  derivatives, 
for  example 

N N 

* 2^F^(1  + In  r^j)  + 4^^  F^(Zj  (A.12) 

1=1  1=1 


(A.6) 

(A.7) 
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they  are  Indeterminate  as  (x^.  z.)  (x.  z. ).  For  removing  this  problem, 

Mo)  I J J i»  1 9 9 

Harder  and  Desmarals'  ' suggested  replacing  the  term  r In  r In  (A. 4) 
with  r In  (r  + e).  Physically  we  apply  distributed  loads  Instead  of  point 
loads  on  the  plate.  The  new  generated  spline  surface  passes  through  N points 
and  has  all  derivatives  everywhere.  The  difficulty  of  using  this  alternative 
Is  to  determine  the  distribution  of  e over  the  spline  surface  and  Harder  and 
Oesmarals^^^^  did  not  give  the  relation  between  e and  other  known  quantities 
explicitly. 

Once  the  values  of  w have  been  found  for  some  constant  x^  (or  z^) 

using  the  surface  spline,  the  derivatives  can  be  obtained  by  using  the  cubic 

linear  spline  because  w Is  now  a function  of  one  variable.  The  cubic  spline. 

In  which  a piecewise,  third-order  polynomial  Is  used  to  describe  a function 

between  two  discrete  points,  has  been  widely  used  to  Interpolate  and  find 

derivatives  of  a function  of  one  variable  because  It  Is  easily  to  calculate  and 

gives  generally  setisfactory  results. 

Once  the  values  of  w and  Its  first  derivatives  dw/dx  and  aw/az  have 

been  found  for  some  constant  x^  (or  z^)  using  the  surface  spline,  the 

second  derivatives  can  be  obtained  by  the  Fourier  series  expansion  method 

(11) 

developed  by  James'  at  Douglas  Aircraft  Company.  The  method  has  the 
advantage  that  It  also  smoothes  the  data  and  the  results  are  not  sensitive  to 
Irregularity  of  the  Input  data. 


APPENDIX  B 


USER'S  MANUAL 

The  present  computer  program  is  written  such  that  various  orthogonal 
coordinate  systems  can  be  used.  In  Case  I,  the  user  provides  the  body  geometry 
and  the  corresponding  inviscid  velocity  components  on  the  body  surface  in  terms 
of  the  Cartesian  coordinate  system.  The  geometrical  parameters  h^,  h2t 
are  computed  internally  for  the  coordinate  system  used  in  this  report.  In 
Case  II,  the  user  supplies  the  geometrical  parameters  and  the  inviscid  velocity 
components  u^  and  along  the  coordinate  lines.  The  input  data  for  each 
case  is  described  below. 

Case  I 

1.  Read  the  values  of  ICASE,  NXT  and  NZT  in  15  Format,  three  variables  on 
one  card. 

ICASE  1 

NXT  Total  number  of  x-stations  for  computing  boundary  layers, 

3 1 NXT  ^ 35 

NZT  Total  number  of  z-stations  for  computing  boundary  layers, 

3 < NZT  i 20 

2.  Read  the  values  of  X(I),  1*1 ,2,. .. ,NXT  in  FIO.O  Format,  one  set  of 
eight  per  card. 

X(I)  Values  of  x where  the  boundary  layers  are  computed  in  meter 

or  ft. 

X(I)  denote  the  Cartesian  coordinates  on  the  body. 

3.  Read  the  values  of  Z(K),  K=1 ,2 NZT  in  FIO.O  Format,  one  set  of 

eight  per  card. 

Z(K)  Values  of  z where  the  boundary  layers  are  computed  In  meter 

or  ft. 

Z(K)  are  in  the  Cartesian  coordinates  on  the  body.  In  order 
to  get  better  results,  a sudden  change  in  z interval  (I.e., 
&Z|^  “ ~ ^K-1  ^ should  be  avoided. 
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4. 


Read  the  values  of  Y,  UE,  VE,  and  WE.  The  variables  are  punched  In 
FIO.O  format,  one  set  of  four  per  card  in  the  following  order: 

{[Y(I.K).UE(I,K),VE(I.K).WE(I,K),K=1,NZT],I-1,NXT} 

Y y-value  at  (x^ , Z|^)  on  the  body  surface. 

bE  Inviscid  velocity  component  along  x-direction  at  (x^,  z^) 

VE  Inviscid  velocity  component  along  y-direction  at  (x^ , Zj^) 

WE  Inviscid  velocity  component  along  z-direction  at  (x^,  Zj^) 

Here  x,  y,  and  z are  in  Cartesian  coordinates,  Y is  in  meter  or 
feet  and  UE,  VE,  and  WE  are  in  m/sec  or  ft/sec. 

5.  Read  TITLE  in  20A4  format.  It  is  used  for  identifying  the  problem  under 
study.  It  should  not  be  greater  than  the  length  of  80  including  blank. 

6.  Read  the  variables,  NXT,  NZT,  NXTRAN,  VI SC,  and  UREF  in  one  card.  The 
first  three  variables  are  in  15  Format,  the  remaining  variables  in  FIO.O 
Format. 

NXT 
NZT 

NXTRAN 

Vise 

UREF 

Case  II 

In  Case  II,  the  input  data  are  almost  the  same  as  for  Case  I except  for 
the  following  changes: 

a.  ICASE  set  equal  to  2. 

b.  X(I)  and  Z(K)  are  given  in  the  corresponding  coordinate  system 
instead  of  the  Cartesian  coordinate  system. 
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Total  number  of  x-stations  for  computing  boundary  layers, 

3 < NXT  < 35 

Total  number  of  z-stations  for  computing  boundary  layers, 

3 NZT  20 

Number  of  x-stations  at  which  the  flow  changes  from  laminar  to 
turbulent.  If  the  flow  is  laminar,  set  the  value  arbitrarily 
greater  than  NXT. 

2 2 

Kinematic  viscosity  of  the  fluid,  in  m /sec  or  ft  /sec. 

A reference  velocity,  in  m/sec  or  ft/sec. 


c.  Replace  the  data  described  In  item  4 by 

{[UE(I,K),WE(I,K),H1(I,K),H2(I,K),CK1(I,K),CK2(I,K),K=1,NZT], 

I»1,NXT) 

The  variables  are  punched  In  FIO.O  format,  six  variables  In  one  card 
UE  Inviscid  velocity  component  alonq  x-directlon 

WE  Inviscid  velocity  component  along  z-directlon 

HI  Metric  coefficient  along  x-direction,  h-| 

H2  Metric  coefficient  along  z-directlon,  hg 

CKl  Geodesic  curvature  coefficient  , see  (2.5) 

CK2  Geodesic  curvature  coefficient  K2,  see  (2.5) 


40 


F 


DISTRIBUTION  LIST 


L , 


\ 

I 


[ 


Technical  Library 
Building  313 

Ballistic  Research  Laboratories 
Aberdeen  Proving  Ground,  MD  21005 

Or.  F.  0.  Bennett 
External  Ballistic  Laboratory 
Ballistic  Research  Laboratories 
Aberdeen  Proving  Ground,  HD  21005 

Mr.  C.  C.  Hudson 
Sandia  Corporation 
Sandia  Base 

Albuquerque,  NM  81115 

Professor  P.  J.  Roache 
Ecodynamics  Research 
Associates,  Inc. 

P.  0.  Box  8172 
Albuquerque,  NM  87108 

Or.  J.  D.  Shreve,  Jr. 

Sandia  Corporation 
Sandia  Base 

Albuquerque,  NM  81115 


Office  of  Naval  Research 
Code  411 

Arlington,  VA  22217 

Office  of  Naval  Research 
Code  421 

Arlington,  VA  22217 

I 

Office  of  Naval  Research 

Cede  438  ^ 

Arlington,  VA  22217  j 

i 

Office  of  Naval  Research  i 

Code  1021 P (ONRL)  | 

Arlington,  VA  22217  (6)  ^ 

Dr.  J.  L.  Potter  I 

Deputy  Director,  Technology  i 

von  Karman  Gas  Dynamics  Facility  | 

Arnold  Air  Force  Station,  TN  37389  | 

Professor  J.  C.  Hu  | 

Georgia  Institute  of  Technology  | 

School  of  Aerospace  Engineering  I 

Atlanta,  GA  30332  I 


Defense  Documentation  Center 
Cameron  Station,  Building  5 
Alexandria,  VA  22314  (12) 

Library 
Naval  Academy 
Annapolis,  iO  21402 

Dr.  G.  H.  Heilmeier 
Director,  Defense  Advanced 
Research  Projects  Agency 
1400  Wilson  Boulevard 
Arlington,  VA  22209 

Mr.  R.  A.  Moore 
Deputy  Director,  Tactical 
Technology  Office 
Defense  Advanced  Research  Projects 
Agency 

1400  Wilson  Boulevard 
Arlington,  VA  22209 


Library  s 

Aerojet-General  Corporation 
6352  North  Irwindale  Avenue 
Azusa,  CA  91702 


NASA  Scientific  and  Technical 
Information  Facility 
P.  0.  Box  8757 

Bal timore/Washington  International 
Airport 

Maryland  21240 

Dr.  S.  A.  Berger 
University  of  California 
Department  of  Mechanical  Engineering 
Berkeley,  CA  94720 

Professor  A.  J.  Chorin 
University  of  California 
Department  of  Mathematics 
Berkeley,  CA  94720 


*1 


A 


J 


Page  2 


Professor  M.  Holt 
University  of  California 
Department  of  Mechanical  Engineering 
Berkeley,  CA  94720 

Dr.  L.  Talbot 
University  of  California 
Department  of  Mechanical  Engineering 
Berkeley.  CA  94720 

Or.  H.  R.  Chaplin 
Code  16 

David  W.  Taylor  Naval  Ship  Research 
and  Development  Center 
Bethesda,  MD  20084 

Code  1800 

David  W.  Taylor  Naval  Ship  Research 
and  Development  Center 
Bethesda , MD  20084 

Code  5643 

David  W.  Taylor  Naval  Ship  Research 
and  Development  Center 
Bethesda , MD  20084 

Dr.  G.  R.  Inger 

Virginia  Polytechnic  Institute 
and  State  University 
Department  of  Aerospace  Engineering 
Blacksburg,  VA  24061 

Professor  A.  H.  Nayfeh 
Virginia  Polytechnic  Institute 
and  State  University 
Department  of  Engineering  Science 
and  Mechanics 
Blacksburg,  VA  24061 

Indiana  University 

School  of  Applied  Mathematics 

Bloomington,  IN  47401 

Oi rector 

Office  of  Naval  Research  Branch 
Office 

495  Summer  Street 
Boston,  MA  02210 

Supervisor,  Technical  Library 
Section 

Thiokol  Chemical  Corporation 
Wasatch  Division 
Brigham  City,  UT  84302 


Dr.  G.  Hall 

State  University  of  New  York  at 
Buffalo 

Faculty  of  Engineering  and 
Applied  Science 

Fluid  and  Thermal  Sciences  Laboratory 
Buffalo,  NY  14214 

Mr.  R.  J.  Vidal 
Cal  span  Corporation 
Aerodynamics  Research  Department 
P.  0.  Box  235 
Buffalo.  NY  14221 

Professor  R.  F.  Probstein 
Massachusetts  Institute  of  Technology 
Department  of  Mechanical  Engineering 
Cambridge,  MA  02139 

Director 

Office  of  Naval  Research  Branch  Office 
536  South  Clark  Street 
Chicago,  IL  60605 

Code  753 

Naval  Weapons  Center 
China  Lake,  CA  93555 

Mr.  J.  Marshall 
Code  4063 

Naval  Weapons  Center 
China  Lake,  CA  93555 

Professor  R.  T.  Davis 
University  of  Cincinnati 
Department  of  Aerospace  Engineering 
and  Applied  Mechanics 
Cincinnati,  OH  45221 

Library  MS  60-3 
NASA  Lewis  Research  Center 
21000  Brookpark  Road 
Cleveland.  OH  44135 

Dr.  J.  0.  Anderson,  Jr. 

Chairman,  Department  of  Aerospace 
Engineering 

College  of  Engineering 
University  of  Maryland 
College  Park,  MO  20742 


-Zt 


Page  3 


Professor  W.  L.  Melnik 
University  of  Maryland 
Department  of  Aerospace  Engineering 
Glenn  L.  Martin  Institute  of 
Technology 

College  Park,  MD  20742 

Professor  0.  Burggraf 
Ohio  State  University 
Department  of  Aeronautical  and 
Astronautical  Engineering 
1314  Kinnear  Road 
Columbus,  OH  43212 

Technical  Library 
Naval  Surface  Weapons  Center 
Dahlgren  Laboratory 
Dahlgren,  VA  22*48 

Dr.  F.  Moore 

Naval  Surface  Weapons  Center 
Dahlgren  Laboratory 
Dahlgren,  VA  22448 

Technical  Library  2-51131 
LTV  Aerospace  Corporation 
P.  0.  Box  5907 
Dallas,  TX  75222 

Library,  United  Aircraft  Corporation 
Research  Laboratories 
Silver  Lane 

East  Hartford,  CT  06108 

Technical  Library 
AVCO-Everett  Research  Laboratory 
2385  Revere  Beach  Parkway 
Everett,  MA  02149 

Professor  6.  Moretti 
Polytechnic  Institute  of  New  York 
Long  Island  Center 
Department  of  Aerospace  Engineering 
and  Applied  Mechanics 
Route  110 

Farmingdale,  NY  11735 

Professor  S.  G.  Rubin 
Polytechnic  Institute  of  New  York 
Long  Island  Center 
Department  of  Aerospace  Engineering 
and  Applied  Mechanics 
Route  110 

Farmingdale,  NY  11735 


Technical  Documents  Center 

Army  Mobility  Equipment  R&D  Center 

Building  315 

Fort  Bel voir,  VA  22060 

Dr.  W.  R.  Briley 

Scientific  Research  Associates,  Inc. 
P.  0.  Box  498 
Glastonbury,  CT  06033 

Library  (MS  185) 

NASA  Langley  Research  Center 
Langley  Station 
Hampton , VA  23665 

Dr.  S.  Nadir 
Northrop  Corporation 
Aircraft  Division 
3901  West  Broadway 
Hawthorne,  CA  90250 

Professor  A.  Chapmann 
Chairman,  Mechanical  Engineering 
Department 

William  M.  Rice  Institute 
Box  1892 

Houston,  TX  77001 

Dr.  F.  Lane 

KLD  Associates,  Inc. 

7 High  Street 
Huntington,  NY  11743 

Technical  Library 
Naval  Ordnance  Station 
Indian  Head,  MD  20640 

Professor  D.  A.  Caughey 
Cornell  University 
Sibley  School  of  Mechanical  and 
Aerospace  Engineering 
Ithaca,  NY  14853 

Professor  E.  L.  Resler 
Cornell  University 
Sibley  School  of  Mechanical  and 
Aerospace  Engineering 
Ithaca,  NY  14853 

Professor  S.  F.  Shen 
Cornell  University 
Sibley  School  of  Mechanical  and 
Aerospace  Engineering 
Ithaca,  NY  14853 


Page  4 


Library 

Midwest  Research  Institute 
425  Volker  Boulevard 
Kansas  City,  MO  64110 

Dr.  M.  M.  Hafez 
Flow  Research,  Inc. 

P.  0.  Box  5040 
Kent,  WA  98031 

Dr.  E.  M.  Murman 
Flow  Research,  Inc. 

P.  0.  Box  5040 
Kent,  WA  98031 

Dr.  S.  A.  Orszag 
Cambridge  Hydrodynamics,  Inc. 

54  Baskin  Road 
Lexington,  MA  02173 

Professor  T.  Cebeci 
California  State  University,  Long 
Beach 

Mechanical  Engineering  Department 
Long  Beach,  CA  90840 

Mr.  J.  L.  Hess 
Douglas  Aircraft  Company 
3855  Lakewood  Boulevard 
Long  Beach,  CA  90808 

Dr.  H.  K.  Cheng 

University  of  Southern  California, 
University  Park 

Department  of  Aerospace  Engineering 
Los  Angeles,  CA  90007 

Professor  J.  D.  Cole 
University  of  California 
Mechanics  and  Structures  Department 
School  of  Engineering  and  Applied 
Science 

Los  Angeles,  CA  90024 

Engineering  Library 
University  of  Southern  California 
Box  77929 

Los  Angeles,  CA  90007 
Dr.  C.  -M.  Ho 

University  of  Southern  California, 
University  Park 

Department  of  Aerospace  Engineering 
Los  Angeles,  CA  90007 


Dr.  T.  D.  Taylor 

The  Aerospace  Corporation 

P.  0.  Box  92957 

Los  Angeles,  CA  90009 

Commandirg  Officer 
Naval  Oroiiance  Station 
Louisville,  KY  40214 

Mr.  B.  H.  Little,  Jr. 
Lockheed-Georgla  Company 
Department  72-74,  Zone  369 
Marietta,  GA  30061 

Or.  C.  Cook 

Stanford  Research  Institute 
Menlo  Park,  CA  94025 

Professor  E.  R.  G.  Eckert 
University  of  Minnesota 
241  Mechanical  Engineering  Building 
Minneapolis,  MN  55455 

Library 

Naval  Postgraduate  School 
Monterey,  CA  93940 

McGill  University 
Supersonic-Gas  Dynamics  Research 
Laboratory 

Department  of  Mechanical  Engineering 
Montreal  12,  Quebec,  Canada 

Librarian 

Engineering  Library,  127-223 
Radio  Corporation  of  America 
Morristown,  NJ  07960 

Dr.  S.  S.  Stahara 

Nielsen  Engineering  & Research,  Inc. 
510  Clyde  Avenue 
Mountain  View,  CA  94043 

Engineering  Societies  Library 
345  East  47th  Street 
New  York,  NY  10017 

Professor  A.  Jameson 
New  York  University 
Courant  Institute  of  Mathematical 
Sciences 

251  Mercer  Street 
New  York,  NY  10012 


I 


Page  5 


Professor  G.  Miller 
New  York  University 
Department  of  Applied  Science 
26-36  Stuyvesant  Street 
New  York,  NY  10003 

Office  of  Naval  Research 
New  York  Area  Office 
715  Broadway  - 5th  Floor 
New  York,  NY  10003 

Dr.  A.  Vagllo-Laurin 
New  York  University 
Department  of  Applied  Science 
26-36  Stuyvesant  Street 
New  York,  NY  10003 

Professor  S.  Weinbaum 
Research  Foundation  of  the  City 
University  of  New  York  on  behalf 
of  the  City  College 
1411  Broadway 
New  York,  NY  10018 

Librarian,  Aeronautical  Library 
National  Research  Council 
Montreal  Road 
Ottawa  7,  Canada 

Lockheed  Missiles  and  Space  Company 
Technical  Information  Center 
3251  Hanover  Street 
Palo  Alto,  CA  94304 

Director 

Office  of  Naval  Research  Branch 
Office 

1030  East  Green  Street 
Pasadena,  CA  91106 

California  Institute  of  Technology 
Engineering  Division 
Pasadena,  CA  91109 

Library 

Jet  Propulsion  Laboratory 
4800  Oak  Grove  Drive 
Pasadena,  CA  91103 

Professor  H.  Llepmann 
California  Institute  of  Technology 
Department  of  Aeronautics 
Pasadena,  CA  91109 


Mr.  L.  I.  Chasen,  MGR-MSD  Lib. 

General  Electric  Company  ] 

Missile  and  Space  Division  i 

P.  0.  Box  8555  ^ 

Philadelphia,  PA  19101  j 

j 

Mr.  P.  Dodge  3 

A1 research  Manufacturing  Company  ] 

of  Arizona 

Division  of  Garrett  Corporation 
402  South  36th  Street 
Phoenix,  AZ  85034 

Technical  Library 
Naval  Missile  Center 
Point  Mugu,  CA  93042 

Professor  S.  Bogdonoff 
Princeton  University 
Gas  Dynamics  Laboratory 
Department  of  Aerospace  and 
Mechanical  Sciences 
Princeton,  NJ  08540 

Professor  S.  I.  Cheng 
Princeton  University 
Department  of  Aerospace  and 
Mechanical  Sciences 
Princeton,  NJ  08540 

Dr.  J.  E.  Yates 

Aeronautical  Research  Associates 
of  Princeton,  Inc. 

50  Washington  Road 
Princeton,  NJ  08540 

Professor  J.  H.  Clarke 
Brown  University 
Division  of  Engineering 
Providence,  RI  02912 

Professor  J.  T.  C.  Liu 
Brown  University 
Division  of  Engineering 
Providence,  RI  02912 

Professor  L.  Sirovich 

Brown  University 

Division  of  Applied  Mathematics 

Providence,  RI  02912 


Dr.  P.  K.  Dai  (Rl/2178) 

TRW  Systems  Group,  Inc. 

One  Space  Park 
Redondo  Beach,  CA  90278 

Redstone  Scientific  Information 
Center 

Chief,  Document  Section 
Army  Missile  Command 
Redstone  Arsenal,  AL  35809 

U.S.  Army  Research  Office 

P.  0.  Box  12211 

Research  Triangle,  NC  27709 

Professor  M.  Lessen 
The  University  of  Rochester 
Department  of  Mechanical  Engineering 
River  Campus  Station 
Rochester,  NY  14627 

Editor,  Applied  Mechanics  Review 
Southwest  Research  Institute 
8500  Culebra  Road 
San  Antonio,  TX  78228 

Library  and  Information  Services 
General  Dynamics-CONVAIR 
P.  0.  Box  1128 
San  Diego,  CA  92112 


Library 

The  Rand  Corporation 
1700  Main  Street 
Santa  Monica,  CA  90401 

Department  Librarian 
University  of  Washington 
Department  of  Aeronautics  and 
Astronautics 
Seattle,  WA  98105 

Dr.  P.  E.  Rubbert 

Boeing  Commercial  Airplane  Company 
P.  0.  Box  3707 
Seattle,  WA  98124 

Mr.  R.  Feldhuhn 
Naval  Surface  Weapons  Center 
White  Oak  Laboratory 
Silver  Spring,  MD  20910 

Dr.  G.  Heiche 

Naval  Surface  Weapons  Center 
Mathematical  Analysis  Branch 
Silver  Spring,  MD  20910 

Librarian 

Naval  Surface  Weapons  Center 
White  Oak  Laboratory 
Silver  Spring,  MD  20910 


Dr.  R.  Magnus 

General  Dynamics-CONVAIR 

Kearny  Mesa  Plant 

P.  0.  Box  80847 

San  Diego,  CA  92138 

Mr.  T.  Brundage 
Defense  Advanced  Research 
Projects  Agency 
Research  and  Development 
Field  Unit 
APO  146,  Box  271 
San  Francisco,  CA  96246 

Office  of  Naval  Research 
San  Francisco  Area  Office 
760  Market  Street  - Room  447 
San  Francisco,  CA  94102 


Dr.  J.  M.  Solomon 
Naval  Surface  Weapons  Center 
White  Oak  Laboratory 
Silver  Spring,  MD  20910 

Professor  J.  H.  Ferziger 
Stanford  University 
Department  of  Mechanical  Engineering 
Stanford , CA  94305 

Professor  K.  Karamcheti 
Stanford  University 
Department  of  Aeronautics  and 
Astronautics 
Stanford,  CA  94305 

Professor  M.  van  Dyke 
Stanford  University 
Department  of  Aeronautics  and 
Astronautics 
Stanford , CA  94305 


:| 


i 


[ 

Page  7 


Engineering  Library 
McDonnell  Douglas  Corporation 
Department  218,  Building  101 
P.  0.  Box  516 
St.  Louis,  MO  63166 

Dr.  R.  J.  Hakkinen 

McDonnell  Douglas  Corporation 

Department  222 

P.  0.  Box  516 

St.  Louis,  MO  63166 

Dr.  R.  P.  Helnisch 
Honeywell,  Inc. 

Systems  and  Research  Division  - 
Aerospace  Defense  Group 
2345  Walnut  Street 
St.  Paul,  MN  55113 

Professor  R.  G.  Stoner 
Arizona  State  University 
Department  of  Physics 
Tempe,  AZ  85721 

Dr.  N,  Malmuth 
Rockwell  International 
Science  Center 
1049  Camino  Dos  Rios 
P.  0.  Box  1085 
Thousand  Oaks,  CA  91360 

Rockwell  International 
Science  Center 
1049  Camino  Dos  Rios 
P.  0.  Box  1085 
Thousand  Oaks,  CA  91360 

The  Library 
University  of  Toronto 
Institute  of  Aerospace  Studies 
Toronto  5,  Canada 

Professor  W.  R.  Sears 
University  of  Arizona 
Aerospace  and  Mechanical  Engineering 
Tucson,  AZ  85721 

Professor  A.  R.  Seebass 
University  of  Arizona 
Department  of  Aerospace  and 
Mechanical  Engineering 
Tucson,  AZ  85721 


Dr.  S.  M.  Yen 
University  of  Illinois 
Coordinated  Science  Laboratory 
Urbana,  IL  61801 

Dr.  K.  T.  Yen 
Code  3015 

Naval  Air  Development  Center 
Warminster,  PA  18974 

Air  Force  Office  of  Scientific 
Research  (SREM) 

Building  1410,  Bolling  AFB 
Washington,  DC  20332 

Chief  of  Research  & Development 
Office  of  Chief  of  Staff 
Department  of  the  Army 
Washington,  DC  20310 

Library  of  Congress 

Science  and  Technology  Division 

Washington,  DC  20540 

Director  of  Research  (Code  RR) 
National  Aeronautics  and 
Space  Administration 
600  Independence  Avenue,  SW 
Washington,  DC  20546 

Library 

National  Bureau  of  Standards 
Washington,  DC  20234 

National  Science  Foundation 
Engineering  Division 
1800  G Street,  NW 
Washington,  DC  20550 

Mr.  W.  Koven  (AIR  03E) 

Naval  Air  Systems  Command 
Washington,  DC  20361 

Mr.  R.  Slewert  (AIR  320D) 

Naval  Air  Systems  Cotimand 
Washington,  DC  20361 

Technical  Library  Division  (AIR  604) 
Naval  Air  Systems  Command 
Washington,  DC  20361 


Page  8 


Code  2627 

Naval  Research  laboratory 
Washington,  DC  20375 

SEA  03512 

Naval  Sea  Systems  Command 
Washington,  DC  20362 

SEA  09G3 

Naval  Sea  Systems  Command 
Washington,  DC  20362 

Dr.  A.  L.  Slafkosky 
Scientific  Advisor 
Commandant  of  the  Marine  Corps 
(Code  AX) 

Washington,  DC  20380 
Director 

Weapons  Systems  Evaluation  Group 
Washington,  DC  20305 

Dr.  P.  Baronti 
General  Applied  Science 
Laboratories,  Inc. 

Merrick  and  Stewart  Avenues 
Westbury,  NY  11590 

Bell  Laboratories 
Whippany  Road 
Whippany,  NJ  07981 

Chief  of  Aerodynamics 
AVCO  Corporation 
Missile  Systems  Division 
201  Lowell  Street 
Wilmington,  MA  01887 

Research  Library 
AVCO  Corporation 
Missile  Systems  Division 
201  Lowell  Street 
Wilmington,  MA  01887 

AFAPL  (APRC) 

AB 

Wright  Patterson,  AFB,  OH  45433 

Dr.  Donald  J.  Harney 
AFFDL/FX 

Wright  Patterson  AFB,  OH  45433 


